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Abstract 

The properties of linear Alfven, slow, and fast magnetoacoustic waves 
for uniform plasmas in relativistic magnetohydrodynamics (MHD) 
are discussed, augmenting the well-known expressions for their phase 
speeds with knowledge on the group speed. A 3 + 1 formalism is 
purposely adopted to make direct comparison with the Newtonian 
MHD limits easier and to stress the graphical representation of their 
anisotropic linear wave properties using the phase and group speed 
diagrams. By drawing these for both the fluid rest frame and for a 
laboratory Lorentzian frame which sees the plasma move with a three- 
velocity having an arbitrary orientation with respect to the magnetic 
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field, a graphical view of the relativistic aberration effects is obtained 
for all three MHD wave families. Moreover, it is confirmed that the 
classical Huygens construction relates the phase and group speed di- 
agram in the usual way, even for the lab frame viewpoint. Since the 
group speed diagrams correspond to exact solutions for initial condi- 
tions corresponding to a localized point perturbation, their formulae 
and geometrical construction can serve to benchmark current high- 
resolution algorithms for numerical relativistic MHD. 
47.35.Tv, 47.75.-|-f, 52.27.Ny, 52.35.Bj 
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I. INTRODUCTION: RELATIVISTIC MHD 



The classical magnetohydrodynamic (MHD) description for the macroscopic 
dynamics of plasmas offers a uniquely powerful, unifying viewpoint on both labora- 
tory and astrophysical plasmas (see e.g. the textbook by Goedbloed and Poedts^). 
The applicability of the MHD description is discussed in many textbooks, and 
it is well known that for most laboratory plasmas, the single fluid ideal or re- 
sistive MHD model eventually needs to be extended towards a multifluid model 
and/or by including important kinetic effects, since its continuum approach to 
plasma modeling neglected e.g. Landau damping as well as many other velocity- 
space dependent physical phenomena. For many astrophysical plasmas, we face 
another shortcoming of the MHD model, namely that it restricts all attention to 
non-relativistic plasma velocities. This latter shortcoming is easily remedied by 
formulating the ideal MHD equations in a frame-invariant relativistic formulation 
within four-dimensional space-time. The mathematical formulation of the rela- 
tivistic ideal MHD equations, and their physical implications, can be found in 
(somewhat more technical) monographs by Lichnerowicz^ and Anile^. The ideal 
relativistic MHD (RMHD) equations combine the full set of Maxwell equations 
with particle and tensorial energy-momentum conservation, under the specific as- 
sumption that the electric field in the co-moving frame vanishes identically, i.e. 
E' = 0, or E = —V X B as expressed in the familiar three-vectors in a fixed lab- 
oratory Lorentzian reference frame. The linear wave dynamics in homogeneous 
relativistic plasmas are also completely known: just as in classical MHD, a ho- 
mogeneous magnetized plasma supports slow and fast magnetoacoustic waves, as 
well as Alfven waves. Expressions for the wave speeds, as obtained for different 
reference frames in relative motion to eachother, have been used as basic ingre- 
dients for the current suite of relativistic magnetofluid codes^"^^. Knowledge of 
the (fastest) characteristic speeds embedded in the relativistic MHD equations is 
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essential to obtain timestep stability limits for explicit time stepping strategies, 
while many modern RMHD codes use more complete information on all linear 
wave dynamics. This information is incorporated in the eigenvectors of the flux 
Jacobian, with the fluxes appearing in the governing hyperbolic equations for the 
conserved variables, typically formulated in a 3 + 1 formulation where a fixed lab- 
oratory frame is selected. In this paper, we (re)collect the basic formulae for linear 
RMHD waves, and present them in a more graphical sense using the classical phase 
and group speed diagram representations, as originally introduced by Friedrichs. 
The relativistic modifications from the classical MHD results are thereby pointed 
out, as well as the MHD wave anisotropics. We give both phase and group speed 
diagrams for all three RMHD wave families, and explain how well-known effects 
for relativistic wave transformations (in particular, relativistic wave aberration) 
connect these diagrams when drawn from differing frames in relative motion. The 
basic connection between phase and group diagrams, by means of a Huygens con- 
struction to obtain group from phase speed diagrams, is demonstrated for cases 
where the point perturbation is seen to move past with arbitrary velocity v. 

Our results complement the wealth of knowledge on linear and nonlinear 
(shock) wave properties for relativistic MHD, as collected in Anile^ and Lich- 
nerowicz^ and in follow-up studies. An example of the latter is the discussion of 
the properties of Alfven waves for both continuous and discontinuous (shock) solu- 
tions by Komissarov^^, emphasizing the wave polarization properties as viewed in 
different Lorentzian reference frames. In a more recent publication^^, a dispersion 
relation for RMHD linear waves in a homogeneous plasma convectcd at uniform 
speed was analyzed in the most tractable special cases (aligned velocity and mag- 
netic field vector). This study adopted a 3-|- 1 formalism, which is the one adhered 
to here as well. It is pointed out how this result fits in the established theoretical 
framework, and how it can then be used to generalize the findings for linear waves 
in relativistic anisotropic MHD formulations^^, with a pressure difference along 
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and perpendicular to the magnetic field. In Sect. II, we start our presentation 
with linear sound wave dynamics in relativistic compressible gases. This is then 
generalized to homogeneous plasmas in Sect. III. 



Since the magnetoacoustic waves are direct generalizations from the gas dy- 
namical sound waves, it is instructive to iterate on the relativistic hydro case first. 
To obtain the relativistic expression for the sound wave speed, one must perform 
the usual linearization of the relativistic hydrodynamic equations about an equi- 
librium configuration. We do this for a homogeneous gas only, adopting a strategy 
which rewrites the governing covariant equations in 4-dimensional space-time to 
equivalent expressions in a 3 + 1 formalism. The 3 + 1 formalism splits temporal 
and spatial derivatives in a fixed Lorentzian reference frame, and shows clearly 
how the equations become more involved with respect to their Galilean versions. 
Also, this 3 + 1 formalism is the one adopted in modern numerical approaches as 
it was already proposed by York^^. When we choose as primitive variables the 
entropy S, the rest frame proper density p, and the velocity v in our laboratory 
frame, and temporarily restrict to a constant polytropic equation of state, where 
S — pp~^, we have as governing primitive variable equations 



II. RELATIVISTIC HYDRO PRELIMINARIES 
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In these equations, c is the light speed, and we introduced 
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where the expressions between brackets denote the Gahlean hmits where the 
Lorentz factor F — > 1 . In this form, no approximations have been made yet (noting 
that the generahty and the elegance of the covariant formulation has been lost), 
and we can clearly identify all terms denoting relativistic corrections with respect 
to the familiar Euler equations. 

The equations (1) can easily be linearized about a static v = 0, uniform gas with 
constant entropy and density S,p. Assuming a plane wave variation exp{—iu!t + 
ik • x) of all linear quantities Si, pi, Vi we arrive at 

LuSi = , 

upl = pk • vi , 

wvi = -^k(57p^-Vi + P^^i) . (3) 

As in the non-relativistic case, this system admits for 5 solutions, where 3 wave 
modes are at marginal frequency cu — 0. These are the entropy wave with arbitrary 
Si but without density or velocity perturbation, together with the two transverse 
translations (shear waves). The physically more interesting modes are compressible 
perturbations, i.e. sound waves, with k- Vi 7^ 0. These obey the dispersion relation 

k'^c^ h ph c^' ^ ' 

While this expression is specific to a case with a polytropic equation of state, it is 
also well-known that for the more appropriate Synge equation of state^^, the sound 
speed Cg is obtained from 

c2 GzG' + G/z ' ^ ' 
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where G{z) — Ks{z) / K2{z) denotes the modified Bessel ratio, and z — pc^/p. The 
derivative in Eq. (5) obeys G' = — 1 — 5G/z, due to recurrence relations. The 
enthalpy is then found from ph — pc^G. The Synge equation of state ensures that 
the effective polytropic index changes smoothly from its classical value 5/3 for a 
monoatomic gas, to the relativistic value 4/3 at ultrarelativistic internal energies. 
In effect, this ensures that the sound speed in relativistic hydro is limited by 
Cg < c/y/S, in the adiabatic case. 

We can also compute the characteristic speeds for relativistic hydro by noting 
that the Eqs. (1) allow to read off the components of the 5x5 coefficient matrix 
W in the quasi- linear form (for ID spatial variation) 

where the primitive variables are collected in V = (S*, p, Vx-, Vy, VzY" ■ When we com- 
pute the 5 eigenvalues A of the W matrix, one obtains the characteristic equation 



(A - 



A '2i\Vrr 9~9" ~\~ 9~9 



= 0. (7) 



Hence, the characteristic speeds either take the value \ — (which obviously 
correspond to the entropy and shear waves from above) , while the sound waves are 
this time found from a seemingly more complicated quadratic expression. Natu- 
rally, both approaches must agree. The key observation is that by computing the 
characteristic speeds from the W matrix, we in fact linearized the equations about 
a moving plasma, and we need to consider how plane waves in the gas rest frame 
transform relativistically to a moving reference frame. This obviously involves a 
Lorentz transformation, which for a frame L' with coordinates (ci', x') moving with 
velocity v with respect to L is given by 



r - 1 

x = X H -— V V • X — Fv i , 



y2 



(8) 



When we momentarily assume that in L', we have a plane wave with variation 
exp^—iuj't' + ik' ■ x'), one finds directly that frame L will still observe a plane wave 
with variation exp(— ia;i + ik • x) with frequency and wavevector given by 



a; = r (u;' + k' • v) , 



k = k' + V 



'oo'T r - 1 

— + k • v) — — 



(9) 



These familiar expressions quantify the relativistic Doppler effect (i.e. the change 
in frequency) and shows that the wave vector is observed to change direction when 
viewed from a moving vantage point. The latter effect is known as relativistic 
wave aberration. The inverse formulae (identical with a;' a;, k' k and a sign 
change for v <-> — v) allow to find the phase speed for frame L given by Vph — u/k 
from the formula 

c2 c2 + r2(^;ph-n-v)'-^;2,- ^ ' 

We introduced the unit vector n — 'k/k. When we invert this formula to find the 
phase speed in frame L, we get the quartic from Eq. (7), in which A and 

<-> n • V. 

To complete this discussion on plane wave propagation in uniform gases for 
special relativity, we now draw the phase and group speed diagrams in different 
reference frames. This is done in Fig. 1, where the z — x plane is drawn. We show 
at top left the phase I'phn/c and group speed Vgr = du/dk in the gas rest frame, 
where we find isotropic sound wave propagation at speed Cg in all directions, and 
where phase and group speed diagrams overlap. We took p — 1 — p and adopted 
a Synge equation of state. The sound wave diagram will always be interior to 
the inner dashed circle which corresponds to the upper limit c/ \/3, while the 
outer circle indicates the light limit. In the right panel, the phase speed diagram 
is plotted as seen in a reference frame where the plane wave emitter is moving 
at velocity v = O.Qce^ along the horizontal 2;-axis. The wave aberration effects 
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deformed the single circle to a kind of double spiral form. In the bottom panels, 
we indicate how a Huygens construction then yields the corresponding group speed 
diagram in the same reference frame. The group diagram is what evolves from a 
point perturbation in a finite time. The final panel just shows this group speed 
diagram in that frame, and demonstrates how the wave front gets "beamed" into an 
asymmetric (about the position of the point source) oval shape. We conveniently 
suppressed the third direction from these diagrams, which would involve revolving 
the shapes about the 2;-axis. To conclude this section on sound waves, we note that 
the group diagram as appropriate for a point perturbation which is seen to pass by 
at velocity v = ve^ can be obtained in various ways. We saw that graphically, it 
is found from a Huygens construction. Secondly, we can note that in frame L', the 
wave fronts are given by the expressions x' — v'^y.{0')t' and z' — v'^ ,^{9')t', where 
the group velocity components are parametrized by the polar angle 9' in this frame. 
When wc transform these expressions to the lab frame L, we find simply that the 
wave fronts are found from 

r (i + ^) 

which is exactly as expected from relativistic velocity addition formulae. Com- 
pletely analogous, we may actually work with the relativistic addition rule gener- 
ally given by 

v,_f(r + (r-i)^) + ^ 



r(i + 



V-V' 

gr 



(12) 



where Vg^ = Cgii'. Using Eq. (9) quantifying the relativistic aberration, we then 
find 



Vgr V 



^9 



1T^) + ^^(1T1) 



(13) 



In this expression, the factor Q is given by 




(14) 



The two sign combinations in Eq. (13) correspond to parametrizing the wave front 
using either the forward, or the backward travehng wave, which is then consecu- 
tively plotted in all 27r directions. The wave front for the group velocity plotted 
in Fig. 1 can be checked to agree with all three (differing) parametrizations: the 
one given in Eq. (11) and both sign combinations of Eq. (13). 



To obtain the propagation speeds for linear waves in relativistic MHD, the 
governing conservation laws in tensorial form can be linearized in space-time. The 
results for the characteristic speeds in RMHD are known and their derivation 
from the covariant equations can e.g. be found in Lichnerowicz^ or Anile"^. Here, 
we revisit these results in a 3 -|- 1 formalism. The algebra involved is then more 
cumbersome (see e.g. the results from Kalra and Gebretsadkan^^), even in the case 
of linearizing about a stationary, homogeneous plasma. This is partly because 
of the wave aberration effects, which we mentioned already for the relativistic 
gas dynamic case, and due to the intrinsic anisotropic nature of the MHD wave 
families. For our purposes, it suffices to note that it is possible to write down 
the equivalent set of equations in a 3 -|- 1 formalism for the primitive variables 
(5, p, V, B), analogous to the Eqs. (1) for gas dynamics. The equation for the 
entropy is identical for ideal RMHD, while the added equations for the magnetic 
field are familiar from the non-relativistic case: 



III. PHASE AND GROUP DIAGRAMS FOR WAVES IN 



HOMOGENEOUS PLASMAS 



f V-B = 0, 



< 




V X (v X B) = . 



(15) 
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The continuity as well as the momentum equation become fairly cumbersome ex- 
pressions in a 3 + 1 split, and we only mention what results from them, after 
linearizing with plane waves exp(— + ik • x). 



A. Rest frame expressions 

The analysis is tractable for the special case of the plasma rest frame, and a very 
elegant means to obtain the relativistic variants of the slow, Alfven and fast wave 
speeds can e.g. be found in the appendix from Komissarov^. For the homogeneous 
plasma rest frame, linearizing and indicating as usual the background (uniform) 
quantities with S, p, B, and the linear variables with Si, pi, Bi, Vi, we get 

ujpi = pk-vi, 

a;Bi = B(k-vi) -vi(k-B), k • Bi = , 

2 2 

uvi = -k (Sjp^-'pi + p^Si) + (k(B ■ Bi) - Bi(k ■ B)) 

+^^^^ {Sip'-'pi + p'Si) . (16) 

We here again temporarily assumed a polytropic equation of state, where the 
specific enthalpy is /i = + jSp'^'^/lj — 1), and introduced the quantity 

w = ph-\ . (17) 

Po 

Note that the relativistic invariant magnetic pressure can be evaluated from the 
rest frame value pmag = B^/2iiq. One can directly compare these expressions (16) 
with the non-relativistic ones, as found in many textbooks, such as Goedbloed 
and Poedts^. Only the momentum equation yields an extra term (the last term is 
purely relativistic, and the coefficients for the other terms are changed to involve 
w). Not surprisingly then, the 7 wave solutions from classical ideal MHD return in 
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slightly modified form. The marginal entropy mode is identical, being the solution 
at a; = for which 5*1 7^ only. The Alfven waves return in a virtually unmodified 
form: they represent solutions with Vi 7^ and Bi 7^ while 

= ^1 = k • vi = k • Bi = B • Bi = B • vi = , (18) 

this time given by the dispersion relation 

They retain their purely field-sampling property familiar from non-relativistic 
MHD and we can express their phase Vph and group velocity Vg^ as follows. For 

that purpose, assuming n = k/Zc and denoting by the angle between n and B, 
we find 

Vph^ScOST?^^ Vgr ^ B 

The phase diagram consists of two circles with diameter B/y/fiQW, and the group 
diagram shows pure pointlike signals along the perturbed fieldline traveling at the 
Alfven speed, as in the classical case. Unlike the sound speed, which is limited to 
be less than c/y/3, the Alfven speed can go up to the light speed (and is equal to 
it for vacuum conditions). 

The compressible modes are obtained from straightforward algebraic manipu- 
lations on Eqs. (16) to the dispersion relation 



4 2 
cu — cu 



(ph^ + — ]+cl^ ^ 



+ kW,^^^0. (21) 



Here we purposely wrote this again in terms of the squared sound speed Cg, while we 
can now also introduce the squared Alfven speed vj^ — B'^c^ /jiQW. The expressions 
are then generally valid, with the expressions for specific enthalpy h and sound 
speed Cg depending on the equation of state. These expressions are also collected 
in several works, e.g. by Komissarov^, in the review on particle acceleration and 
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relativistic shocks by Kirk and Duffy^^, and in the early work by Majorana and 
Anile^°. As a result, the phase speeds for slow and fast wave families are found 
from 

Vph 



Vph 

c c 



>^^« + 4|Vl + 5cos2^±a. 

rt^ j 



\^ 2 \ w 

Here, the symbols 5 and a express the following dimensionless ratios 



(22) 



The symbol a follows from 



(23) 



= (l + 5 cos^ -df -(J cos^ ■& . 



(24) 



Noting that ph/w — 1 — v\/c^, the phase speed for purely parallel propagation 
reduces to the same expression found in non- relativistic MHD, where we have 



1 



1 ± 



1 - 



(25) 



The group speed is then written in terms of the orthogonal directions n = k/A; 
and t = [(B/i?) x n] x n as 



Ygr 

c 



c 



nit 



[a ^25 (a± (l + ^cos^ ^9))] sin^9cos^9 



(26) 



2(l + 5cos2'(?±a)a 

These can directly be compared with the non-relativistic expressions (see again 
e.g. Goedbloed and Pocdts^) and we can note that all relativistic effects are due 
to the parameter 5, together with the fact that both sound and Alfven speeds get 
relativistic corrections. As representative examples, we show in Figs. 2-4 the phase 
and group diagrams appropriate for the plasma rest frame, for three cases with 
(again with units making c—1) 

• p=l, p=l, B = l for which Cg = 0.558 and va = 0.431. This is 
shown in Fig. 2, and corresponds to a relativistically hot gas with the sound 
speed approaching its maximal value. The Friedrichs diagrams resemble the 
classical Newtonian versions. 



13 



• p — 1, p — 1, and B — 3 ior which Cg — 0.558 and va — 0.820. This is shown 
in Fig. 3, and corresponds to a hot gas with a very strong magnetic field. 

• p = 0.01, p = 0.001, and 5 = 1 for which Cg = 0.349 and va = 0.994. 
This case has the Alfven speed nearly reaching the light limit and is for a 
cold plasma. Note how both Alfven and slow group diagrams correspond to 
extreme anisotropic, field-sampling perturbations. 

All these assume a Synge equation of state. We can note that the phase speeds obey 
a strict ordering where slow (S), Alfven (A) and fast (F) obey < t'ph,s < Vph.A < 
Vph,F < c for any direction n, and in the plasma rest frame appear symmetric 
about the zero value for forward and backward propagating waves. The fact that 
these rest frame Friedrichs diagrams generalize their Newtonian counterparts was 
already pointed out in the early lecture notes by Newcomb^^. 

B. Lab frame expressions 

To make the RMHD wave dynamics complete, we still need to show phase and 
group diagrams for inertial frames in relative motion to the homogeneous plasma 
rest frame. This can indeed be done by linearizing about a moving (i.e. stationary) 
uniform plasma, as was pursued in Kalra and Gebretsadkan^"^ for the special case 
where the movement v was aligned with the uniform field B, but this is by far 
the most algebraically complex means to do so. It must be realized that the more 
general result (for arbitrary orientation between v and B) was already given in 
the textbook by Lichnerowicz^ and Anile^, obtained from the more appropriate 
covariant form (i.e. valid in any reference frame). Still, because much intuition 
has been built up in classical MHD in a 3 + 1 framework, and modern special 
relativistic codes exploit this as well, we now discuss the RMHD waves in the lab 
frame using the 3 + 1 split. In fact, we then just need one additional equation. 
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relating the three- vector magnetic fields in L' (moving away with speed v) and L, 
by the well known Lorentz formula (here for ideal RMHD) 

. B vB r 

If we indicate the rest frame as L', we know from Eq. (19) that (fph/c)^ = 
(B' • n')^ ///o^y'■ Then, combining Eq. (27), the wavevector transformation from 
Eq. (9) and the phase speed relation in Eq. (10), we obtain the phase speed for 
Alfven waves in a frame where the plane wave emitters move away with velocity 
V as 

\ 



Vph _ t^ph^ _ " • ^ _^ 

C C \ C ]72 



vB 



n. (28) 



The magnetic pressure can be obtained in the lab frame from 



2P„..= ^ + ^. (29) 

It is straightforward to show that Eq. (28) gives once more two circles through 
the origin, and that correspondingly, the group speed diagram still consists of 
two point perturbations, this time seen to propagate on the fieldline which was 
originally through the origin at the time of the initial pointlike perturbation. The 
group speed in the lab frame can be obtained by using the rest frame expression 
Vgj./c = ±B'/ \J \Xqw' and use velocity addition rule in Eq. (12) to get 

^ = ^T ° (30) 

The phase and group speed diagrams for the Alfven waves, in a case where B 
is oriented along the z-axis and the source moves away with a velocity v = 
0.9[sin(7r/4)e2. + cos(7r/4)e2] (without loss of generality taken in the x — z plane), 
are shown in Fig. 5, for similar values as used in Figs. 2-4. The diameter of the 
circles (actually spheres, but we suppressed the third dimension in the plots) that 
define the phase diagram is set by the forward and backward group velocity, and 
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clearly, the part of v perpendicular to B merely displaces the field fine originally 
through the origin (in accord with the frozen-in nature of ideal MHD), on which 
the pointlike Alfven signals propagate. 

Finally, for the slow and fast wave families, the phase speed in the lab frame 
is again best computed from transforming the rest frame expression in Eq. (21) to 
the lab frame. This was already pointed out by Komissarov^, and in doing so we 
obtain the quartic 

,k - 4) - - (i - f ) {r= (,,4 + 2,^^) (^ - 



-4 



Y (- ^ 1 _ ^LJC] _ " ^ 
\c ^/JI^) \ c c ) T^/JI^ 



= 0. (31) 

This expression has been known for a long time^'^° (but in perhaps less transparant 
notation than the one adhered to here), and is the one used in RMHD codes to 
compute the fastest wave speed from, needed for explicit time integration strate- 
gies^^^'^^. The cold plasma limit is obtained for Cg — 0, and has degenerate slow 
waves. 

Wc plot in Fig. 6 the phase diagrams for the four solutions from the quartic 
Eq. (31), as they are seen in the lab frame, for the two cases from Figs. 3-4. The 
Alfven phase speeds are shown as well. We plotted the phase speeds only in the 
x — z plane (plotted is thus Vph = VphU by varying n'va.x — z over 27r), and it should 
be noted that the suppressed third dimension is no longer a mere revolution about 
some axis as soon as the velocity v is not ahgned with B. For the particular case 
where v is aligned with B (the case in the study by Kalra and Gebretsadkan^"^), it is 
merely to be rotated about the 2;-axis. In any case, the visual representation of the 
phase diagram is arguably the most intuitive way to illustrate the combined effects 
of the inherent anisotropic nature of MHD wave propagation (even in uniform 
media), with the added complexity brought in by relativistic wave aberration. As 
in the case of the sound waves, as we vary the unit vector n over 27r (or the full 
47r steradian), the pair of forward/backward slow waves will trace out the same 

16 



curve (surface) in a different parametrization, and similarly for the fast wave pair. 
It is seen from Fig. 6 that the fast wave phase diagram deforms much like the 
gas dynamic sound wave. As the transition from the rest frame to the lab frame 
involves nothing else than what is given by Lorentz transformations, the ordering 
of the forward and backward waves remains in any direction, but the symmetry 
about the zero velocity is obviously lost, and we get 

-C < tiph.F- < ^'ph,A- < ^'ph,S- < Vph,S+ < Vph,A+ < Vph,F+ < C . (32) 

The two examples drawn in Fig. 6 demonstrate that the manner in which the 
Lorentz transformation deforms the rest frame version is truly nontrivial, and can 
give rise to easthetically appealing diagrams. 

Finally, we demonstrate in Fig. 7 how one can again use a mere Huygens con- 
struction to obtain the group speed diagram appropriate for a point perturbation 
that is seen to pass by at velocity v. The left panel only shows the leftward part 
of the two cusp-shaped slow signals that arise from performing the Huygens con- 
struction on the slow phase curve drawn in the left panel of Fig. 6. The right panel 
repeats the fast phase diagram, and demonstrates the group diagram which is now 
a fully three-dimensional surface (in the figure seen in the x — z cross-section) . In 
the special case with aligned v and B, we also conclude that the transformation 
from Eq. (11) can be used directly on the group diagram found in the rest frame. 
The result of this for the cases from Fig. 2-4 is shown in Fig. 8, when the velocity 
is at 90% of the hght speed. Closed form expressions for the group speed diagram 
in the lab frame arc not given explicitly for arbitrary orientation between v and 
B, but the procedure to obtain them is now clear in principle (but an algebraic 
nightmare): we need to perform relativistic speed addition using the rest frame 
expressions from Eq. (26), and thereby transform not just the normal, but also 
the tangential components of the magnetic field to the wave front from frame V 
to the lab frame L. The geometric Huygens construction is a much easier means 
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to realize this. 

IV. CONCLUSIONS 

We provided graphical insight into the phase and group speed properties for 
RMHD waves in uniform, stationary plasmas. The added value of knowing the 
group speed diagrams in different laboratory frames is foremost important to test 
the current suites of RMHD codes: we here pointed out how to obtain the exact 
solution (for Alfven in closed form, and for fast and slow, graphically) to the initial 
value problem where a pointlike (linear) perturbation is set up in an overall station- 
ary, homogeneous plasma. Varying the governing equilibrium parameters (plasma 
beta, looking at the special case where Cg — va, from warm to cold plasma condi- 
tions, from low to high advection speeds), one can found out more about intrinsic 
limitations in the numerical treatment (numerical diffusion, dispersion, stability, 
etc.). The fact that pointlike perturbations in stationary plasmas will give group 
diagrams that can be very demanding on spatial resolution (due to the wave aber- 
ration effects combined with MHD wave anisotropy) , is yet again a clear indication 
that RMHD codes will need some form of grid-adaptivity, to handle complex 3D 
RMHD astrophysical problems. The linear wave properties discussed here are also 
relevant to better appreciate the knowledge of discontinuous (shock) wave solutions 
allowed by the RMHD equations^^'^'^'^°'^^'^^, or the hnear wave modifications en- 
countered in fluid models for relativistic plasmas that invoke anisotropic pressure^^ 
or do not have sufficient coUisionality to justify a RMHD viewpoint^^'^^. Finally, 
knowledge of the wave properties in uniform media is indispensible to appreciate 
the signiflcant modifications encountered when diagnosing waves and instabilities 
in non-uniform, relativistic MHD equilibrium configurations. This would lead to 
MHD spectroscopy^^ for accretion disks about compact objects, or for relativistic 
jets. 
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Figure Captions: 

Figure 1: Group and phase diagrams for a relativistic hydro case with p = 
1 = p and a Synge equation of state. Top left: in the gas rest frame, where the 
circular group and phase diagram coincide, and lie within the limit 1 / y/Sc (dotted 
curve). The light limit is indicated as well (taking c = 1). Top right: the phase 
diagram from a lab viewpoint which saw the source pass at v = O.Qe^. Bottom left: 
a Huygens construction in the same lab viewpoint, to obtain the group diagram 
for the lab frame. The source position is indicated by a small square. Bottom 
right: the group diagram for the lab frame, giving the wave fronts as they appear 
from a point perturbation passing by at v = O.Qe^. 

Figure 2: Phase (left) and group (right) diagram in the plasma rest frame, 
for a case with p = 1 = p = B and a Synge EOS. In the left panel, the maximal 
sound speed limit is indicated as well (dotted curve), and both panels also show 
the light limit. 

Figure 3: Phase (left) and group (right) diagram in the plasma rest frame, for 
a case with p = 1 = p and B = 3 and a Synge EOS. In the left panel, the maximal 
sound speed limit is indicated as well (dotted curve), and both panels also show 
the light limit. 

Figure 4: Phase (left) and group (right) diagram in the plasma rest frame, 
for a case with p = 0.01, p = 0.001 and B = 1 and a Synge EOS. In the left panel, 
the maximal sound speed limit is indicated as well (dotted curve) , and both panels 
also show the light hmit. 

Figure 5: Phase (left) and group (right) diagram for the Alfven waves only, 
in the lab frame where the perturbation of the central (horizontal) field line is 
seen to move at v = 0.9 [sin(7r/4)ej; + cos(7r/4)e2;]. Top panel is for a case where 
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p— 1— p — B{asin Fig. 2), bottom panel for a case where p — 0.01, p — 0.001 
and B = 1 (as in Fig. 4). 

Figure 6: Phase diagram for fast, Alfven, and slow waves, in the lab frame 
where the perturbation of the central (horizontal) field line is seen to move at 
V = 0.9 [sin(7r/4)e^ + cos(7r/4)e2]. The case at left is p = 1 = p and B — 3 {as in 
Fig. 3), at right we have p = 0.01, p = 0.001 and B = 1 {as in Fig. 4). 

Figure 7: Group diagram for slow (left) and fast (right) waves, in the lab 
frame as obtained by means of a Huygens construction. The perturbation of the 

central (horizontal) field line is seen to move at v = 0.9 [sin(7r/4)e2; + cos(7r/4)e;5]. 
The case is p = 1 = p and B — 3 {as in Fig. 3 and left panel of Fig. 6). 

Figure 8: Group diagram for all waves in the lab frame with perturbation 
seen to move at v = O.Qe^. The case is at left p — 1 — p — B (as in Fig. 2) and at 
right p = 0.01, p = 0.001 and £ = 1 (as in Fig. 4). 
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